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We perform molecular dynamics simulations using the extended simple point charge 
SPC/E water model in order to investigate the dynamical behavior of supercooled-stretched 
water. We focus on the behavior of the translational diffusion coefficient, the orientational 
relaxation time, and the local hydrogen bond network. Decreasing density or pressure 
along an isothermal path, there is a mininum in the diffusion coefficient and a maximum 
in the orientational relaxation time, suggesting an initial enhancement and subsequent 
breakdown of the tetrahedral structure and of the hydrogen bond network as the density 
decreases. The analysis of the tetrahedral structure of the nearest neighbors help to clar- 
ify the relationship between the local structural changes and the system dynamics. We 
also find that the product of diffusion coefficient and relaxation time is nearly constant, 
indicating that the anomalous behavior observed in the rotational and translational dif- 
fusion arise from the same microscopic mechanism. 



1 Introduction 

Water is one of the most intriguing subjects of research from the physical, chemical or 
biological point of view. Due to the ubiquity and anomalous thermodynamic behavior of 
water, it has been studied for centuries, and even now unexpected phenomena continue 
to be discovered [^. 

Ordinary liquid water is highly locally ordered, with each molecule having almost ex- 
actly four tetrahedrally coordinated neighbors forming a connected random network of 
hydrogen bonds well above the percolation threshold. Despite possible expectations that 
water might behave like a stable network, the hydrogen bond network restructures itself 
in picosecond time scale. This indicates that the random tetrahedral network is not per- 
fect, but contains some sort of structural defects, what it is supported by experiments. 
The local structure is sensitive to variable temperature or pressure, which can result in 



modifications such as distortion of tlie liydrogen bond angles and distances, as well as 
increase (or decrease) in the number of neighbors. These changes in local order have been 
shown to be intimately connected with the thermodynamic anomalies 0. This is the 
case of the existence of a maximum in the density at 4°C, unusually high melting, boiling 
and critical points, existence of a minimum in the constant pressure heat capacity, among 
several others. 



These local network defects are also important for the mobility of water. In contrast to 
conventional wisdom, the increase of local density (and therefore the increase in pres- 
sure) leads to increased mobility in ambient water, a result of weakening and possible 
destructions of hydrogen bonds that favor an open structure, as shown by experimental 
observations [0, Q and simulations 0, ^, 01- At sufficiently large pressure, steric constraints 
become significant and mobility decreases on increasing density, as in common liquids. As 
a result, along each isotherm, there is a maximum of the diffusion coefficient as a function 
of density. 

If the pressure is decreased it appears an opposite behavior, and if we analyze the diffusion 
coefficient along isotherms, a minimum of the diffusion coefficient can be found at a given 
low density [^, |^ As a consequence, a minimum in the rotational correlation time is also 
expected. 



The behavior of water over a wide range of pressures and temperatures can be summa- 
rized in a coherent picture 0. The region in the phase diagram in which thermodynamic 
anomalies occur is entirely inside the region of kinetic or dynamical anomalies which in 
its turn is contained in the region of structural anomalies. This picture suggests that all 
the anomalous behavior of water can be explained on basis of its structure. In this case, 
within the region of structural anomalies the orientational and translational order should 
be inter-dependent In order to check this hypothesis, here we study the behavior of 
the orientational relaxation and diffusion coefficient of stretched supercooled water. By 
comparing these two quantities, as well as by studying the local structure of water [lU 



such as revealed by the angular distribution of neighbors and hydrogen bonds, we will be 
able to understand how the structure affects the mobility. 



2 Methods 



We performed molecular dynamics simulations using 216 water molecules described by 
the extended simple point charge (SPC/E) model |Tl]], in the canonical ensemble (NVT), 
in a cubic simulation box using periodic boundary conditions. The Berendsen method 
for rescaling the velocities was applied |T^, the electrostatic interactions were calculated 
using the reaction field method |jl3|] with cutoff of 0.79 nm. The equations of motion were 
solved using the SHAKE algorithm |15] with time steps of 1.0 fs for T > 210 K and 
2.0 fs for T = 210 K. 



The diffusion coefficient D was calculated from the asymptotic slope of the time depen- 
dence of the mean square displacement. The orientational relaxation was analyzed using 
the rotational autocorrelation functions IM: 



e) = (P, [e(t) ■ e(0)]) 



(1) 



where e is a chosen unity vector describing the orientation of the molecule and Pi is the 
i-th order Legendre Polynomial. We restrict ourselves to the analysis of the first two 
Legendre Polynomials, using two choices of reference vectors: (1) the dipole vector and 
(2) a vector perpendicular to the plane of the molecule. The correlation functions were 



fitted to a biexponential decay function |16|: 



C = aoexp{-bt'^ /2) + aiexp^-t/r^) + anexp{-t/T^^) (2) 

but in the most of the cases taking only one exponential yielded a reasonable result, and 
thus only one relaxation time was determined. 

In order to understand the effect of the structure on the mobility we also carried out a 
detailed analysis of the local structure of water |TD[: the distribution of the hydrogen 



bond angle (0-H ■ ■ ■ O) and the distribution of the angle between water molecules linked 
to a central water molecule (that is, the angle O ■ ■ ■ O ■ ■ ■ O). 



3 Results 

Our results are summarized in Table 1. The presence of a second relaxation time (the 
faster one) is found mainly (but not only) at the higher densities, in agreement with the 
interpretation that these modes are linked to a disturbance of a tetrahedral local struc- 



ture [|16l. Since the contribution of this faster relaxation time is almost negligible, we 



show only the slow relaxation times in the Table 1. 



Table 1: Results of the simulations. Tx means the orientational relaxation time obtained 
by exponential fit of the y-th order Legendre Polynomial correlation using the x-th vector 
(1 = dipole, 2 = perpendicular) 
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The diffusion coefficient, as illustrated in Figure 1, has a minimum at low densities If 
we take into account the results at high pressures 0, a sigmoidal behavior is clear. The 
lines of minima and maxima of the diffusion coefficient are the borderlines of the region 
of kinetic or dynamic anomalies, and it can be shown that the thermodynamic anomalies 
occur inside the region of dynamic anomalies [0, |]. 

The minimum in the diffusion coefficient is located at a density similar to ice Ih- It is 
interesting also to analyze, whether the orientational relaxation has a similar dependence 
on density. 
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Fig 1. Diffusion Coefficient versus Density along 
isotherms: Detail of the minimum, from top to bottom 
T = 250 K, 240 K and T = 230 K 



The analysis of the reorientational correlation times show peculiar aspects. First, we can 
look (see Table 1) at the values of the relaxation times for each chosen vector in order to 
see some trends. At high temperatures, the relaxation times using the first order Legendre 
polynomial has nearly three times the value of the corresponding second order Legendre 
polynomial. This ratio decays to about two or even less at lower temperatures. This, 
however, can be attributed to enhancement of the fluctuations at the supercooled region. 



According to the Debye rotational diffusion model, this ratio should be indeed three [|17 



The agreement between our results and the theory is quite striking since for the mobility 
of water in its own medium, the hydrodynamic description should not be valid. 



The behavior of the relaxation times with the density, as can be seen in Table 1, was 
also analyzed for the several vectors using either first or second order Legendre polyno- 
mials. At high temperatures, T = 280 K and T = 260 K, the relaxations times show 
no specific trend with the density. There is a small decrease at low densities, but this 
decrease can be attributed to statistical fluctuations. At lower temperatures a maximum 



is observed at about the same density as the translational diffusion coefficients display 
their minima p, 0. This result shows that the enhancement of the water structure has 
both influence in the translational as well as in the rotational diffusion, as expected. 
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Fig 2. Relaxation times obtained from orientational cor- 
relation as discussed in the text, for simulations at T = 
240 K and several densities 

For studying the mobility of large molecules in a solvent made of small particles, hy- 
drodynamic arguments are valid. In that case, the product of diffusion coefficient and 
orientational relaxation time is a constant, independent of density, viscosity and temper- 
ature . For studying the mobility of water in its own solution, neither the continuous 
description of the solvent nor modeling the water molecule as compact object are appro- 
priated, so the hydrodynamic description is expected not to be valid. Nevertheless, the 
similarity between the results for diffusion coefficient and relaxation times seems to indi- 
cate that these two quantities are related. Therefore we computed the product of -D x r. 
For all simulations (see Figure 3) this product remains almost constant what suggests 
that the translational and rotational mobility are controlled by a common mechanism. 
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Fig. 3 The product of diffusion coefficient and relaxation 
time is nearly constant, independent of the temperature 
and density 

So if not due to hydrodynamics, how can one explain that the correlations between D 
and r? Both anomalies in the translational and rotational diffusion are located inside 
the region of structural anomalies what seems to indicate that the structure controls 
the mobility. Even though we have no definitive proof of this statement, we have some 
clues. A detailed analysis of the local structure of water reveals the enhancement of the 
tetrahedrality at low temperatures and ice-like densities. This can be detected from the 
distribution of the hydrogen bond angle and the distribution of the angle between neigh- 
bors O ■ ■ ■ O ■ ■ ■ O illustrated as a function of density for T = 240 K in Figures 4a,4b and 
as a function of temperature for density 0.925 g cm~'^ in Figures 5a,5b. 

The first clue comes from the O ■ ■ ■ O ■ ■ ■ O angle distribution. At low densities, it is 
centered around the tetrahedral angle, (109.5") confirming the enhancement of the ice-like 
structure. However, a small peak at angle 60 ° is found. This peak increases slightly with 
increasing densities until the density (pm = 1.125 g/cm^) and it is related to the presence 
of a fifth neighbor molecule in the first coordination shell. At higher densities, a shoul- 
der also appears around 90 °, that may be related to the appearance of a 6th neighbor. 
Therefore, at very high densities, the O • ■ • O ■ ■ ■ O distribution looks very different with 
strong distortions in the structure. 

The second clue, comes from the hydrogen bond angle distribution. The effect of the 
density is negligible: at low densities all curves almost collapse in a single curve with only 
one peak at low angles, indicating a rigid structure. Only at high densities the distortions 
become strong. 




angle (degrees) angle (degrees) 



Fig. 4 (a) Hydrogen bond angle distribution and (b) neighbors O ■ ■ ■ O ■ ■ ■ O angle 
distribution for several densities at the temperature T = 240 K. Shown are the results 
for p = 0.875 g (Q), 0.900 (□), 0.925 (o), 1.000 (^) and 1.125 g cm''^ (A) 
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Fig. 5 The same as Fig. 4, but for several temperatures at fixed p = 0.925 g cm ^. 
Shown are the results for T = 280 K (0), T = 240 K (□) and T = 210 K (^) 

As shown in Figure 5a, the increase of temperature leads to a broader hydrogen bond dis- 
tribution, but peaked in the same value, indicating almost no distortion, apart from the 
thermal fluctuation. Figure 5b illustrates the O ■ ■ ■ O ■ ■ ■ O angle distribution. One can 
observe that the presence of interstitial water is not strongly affected by the temperature. 
The increase of temperature has, as only effect, a thermal broadening of the distribution. 



4 Conclusions 



We have studied the mobihty of SPC/E water. In contrast with the behavior of a normal 
fluid, the diffusion coefficient of liquid water increases with the density in a broad range 
of thermodynamic conditions, and exhibits a minimum at the ice-like density and a maxi- 
mum at pm- In compass, the rotational diffusion times increase as the density is decreased 
and has a maximum at the ice-like density. Therefore, in this supercooled stretched re- 
gion, the translational and rotational mobilities are not independent, but correlated. The 
product of D X r is a constant in the region of supercooled stretched water. 

Since the hydrodynamic arguments are not valid for describing the mobility of water in 
its own medium, we look for another mechanism for this correlation between D and r. 
The analysis of the O ■ ■ ■ O • ■ ■ O and hydrogen bond angle distributions in the stretched 
region indicates that the increase in pressure disturbs the structure by inclusion of inter- 
stitial molecules but does not cause significative distortions in the H-bonds. This shows 
that the increase of the number of water molecules found in the first coordination shell 
and the fraction of water molecules with more than four intact H-bonds correlates with 
the increase of the mobility. This suggests that the structural defects are directly associ- 
ated with the anomalous behavior of translational and rotational diffusion. 
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